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SUMMARY 

SPUROPT,  a  computer  program  developed  for  optimizing  spur  gear  designs,  was  updated 
by  installing  a  new  subroutine  that  uses  AGMA  908- B89  standards  to  calculate  the  J-factor  for 
determining  tooth-bending  stress.  The  updated  SPUROPT  program  was  then  used  to  optimize 
a  spur  gear  set  for  either  maximum  fatigue  life,  minimum  dynamic  load,  or  minimum  weight. 

All  calculations  were  made  with  constraints  on  as  many  as  13  parameters  by  using  three  design 
variables:  the  number  of  teeth,  diametral  pitch,  and  tooth-face  width.  Results  depended  largelv 
on  constraint  values.  When  ihe  limiting  bending  stress  was  set  at  a  high  value,  the  optimal 
solution  was  the  highest  allowable  number  of  teeth.  When  the  allowable  bending  stress  was 
lowered,  the  optimal  solution  moved  toward  the  fewest  number  of  teeth  permitted.  Final  results 
were  also  affected  by  the  amount  of  transmission  error.  A  lower  error  permitted  a  higher  num¬ 
ber  of  teeth. 


INTRODUCTION 

Gears  are  an  important  component  of  many  power  transmission  systems.  Their  design  can 
greatly  influence  durability,  reliability,  and  maintainability  of  r  system,  aK  well  as  its  overall 
weight.  As  a  result,  it  is  becoming  increasingly  important  to  develop  the  best  gear  design  for  a 
given  situation,  as  system  design  criteria  become  increasingly  stringent.  One  technique  for  pro¬ 
viding  reliable  design  methodology  is  to  use  numerical  optimization  methods  along  with  an  ana¬ 
lytical  gear  code.  A  study  by  Vanderplaats  resulted  in  a  spur  gear  design  computer  program 
called  SPUROPT.  This  program,  fully  described  in  reference  1,  combined  spur  gear  analysis 
routines  with  the  COPES/ADS  optimization  program  (ref.  2).  In  reference  1  the  example 
problem  for  maximizing  surface  fatigue  life  indicated  an  improvement  of  almost  200  percent  for 
the  optimized  gear  design  over  the  original  design.  This  was  a  single  calculation,  however.  For 
a  more  complete  evaluation  of  SPUROPT,  the  program  needed  to  be  exercised  over  a  range  of 
variables.  This  report  evaluates  the  results  of  running  the  computer  program  SPUROPT  over  a 
range  of  variables. 

PROGRAM  SPUROPT 

SPUROPT,  described  more  fully  in  reference  1,  is  composed  of  two  parts.  One  part  con¬ 
sists  of  analytical  routines  that  calculate  such  properties  as  tooth-bending  stress,  surface  stress, 
dynamic  loading,  scoring  criteria,  flash  temperature,  weight,  and  fatigue  life.  The  second  part 
contains  the  COPES/ADS  optimizing  routines  (refs.  1  and  2).  This  part  of  the  program  seeks  to 
determine  the  proper  combination  of  values  of  the  design  variables  (such  as  number  of  teeth  and 
face  width)  that  results  in  maximization  (or  minimization)  of  an  objective  function  (e.g.,  fatigue 
life)  while  other  properties  such  as  tooth-bending  stress  or  gear  weight  are  subject  to  constraints. 
For  further  discussions  about  numerical  optimization  and  the  structure  of  the  SPUROPT 
program,  see  r?ff‘r*‘r>ro  1 


The  example  problem  in  reference  1  maximized  the  fatigue  life  of  a  set  of  spur  gears  for 
specific  material  and  lubrication  constants,  and  for  a  given  operating  condition  (load,  speed,  or 
temperature),  by  constraining  design  variables  such  as  number  of  teeth,  module  (diametral 
pitch),  and  tooth-face  width.  The  program  was  subject  to  limitations  on  bending  stress,  Hertz 
or  contact  stress,  flash  temperature,  contact  ratio,  center  distance,  gear  weight,  and  tooth-width- 
to-diameter  ratio.  This  example  problem  was  the  starting  point  for  all  computer  runs  described 
in  this  paper. 


RESULTS  AND  DISCUSSION 

Initially,  we  ran  the  original  SPUROPT  program  for  the  same  spur  gear  life-maximization 
problem  as  given  in  reference  1.  Typical  gear  input  data  are  shown  in  table  I,  and  a  sample 
input  data  set  is  shown  in  table  II.  A  complete  explanation  of  required  input  data  is  given  in 
references  1  and  2.  The  program  requires  engineering  units  for  input  values  and  provides  output 
values  in  engineering  units.  All  SI  values  shown  in  this  report  were  obtained  by  conversion. 
During  these  initial  runs,  some  parameters  were  changed  (e.g.,  shaft  diameter  and  bending-stress 
limits)  to  determine  the  effects  of  these  variables.  Then,  the  program  was  updated  by  substitut¬ 
ing  a  different  routine  for  the  J-factor  calculation  and  was  run  again  for  evaluation.  This 
activity  is  discussed  in  the  following  sections. 


Initial  Runs 

The  results  of  the  initial  runs  showed  that  gear  fatigue  life  for  the  final  optimal  values  of 
the  design  variables  increased  over  that  life  obtained  with  the  original  design.  In  watching  the 
optimization  process,  we  observed  interesting  intermediate  results  from  SPUROPT.  Figure  1 
shows  the  values  of  the  design  variables,  as  controlled  by  the  optimizer,  and  values  of  the  objec¬ 
tive  function  (gear  fatigue  life)  for  each  trial  calculation  during  solution  of  the  life-maximization 
problem.  After  the  30th  trial,  there  is  a  large  change  in  tut  number  of  i,eetu  (fig.  1(a))  as  well 
as  a  steady  improvement  in  fatigue  life  (fig.  1(d)). 

Summary  results  for  a  typical  run  are  given  in  table  III,  which  also  includes  results  for  the 
initial  design  values  (table  I)  as  well  as  bounds  placed  on  the  problem.  A  study  of  the  output 
data  shows  very  small  values  for  the  tooth-bending  stress  for  some  trials.  This  low  bending 
stress  was  caused  by  some  unusual  values  for  the  J-factor.  Figure  2  shows  values  of  the  J-factor 
as  calculated  by  the  program  for  each  trial.  A  comparison  of  figure  2  with  figure  1(a)  shows 
that  the  unusual  values  occurred  with  the  higher  numbers  of  teeth.  Note  that  the  J-factor 
should  always  be  positive  and  less  than  1. 

J-factor  calculation  in  SPUROPT  was  based  on  a  model  presented  in  reference  3.  Similar 
unusual  values  of  the  J-factor  have  been  noted  previously  (Errichello.  R.,  1989,  Geartech, 
Albany,  CA,  letter  to  Dennis  Townsend,  NASA  Lewis  Research  Center)  for  an  equivalent  set  of 
equations.  To  obtain  better  J-factor  values,  wc  decided  to  replace  the  subroutine  GEOFAC 
(ref.  1),  which  was  used  to  calculate  the  J-factor,  with  equations  shown  in  AGMA  908-B89 
(ref.  4).  Thus,  a  new  subroutine  called  JFACTR  was  written  and  installed.  Subsequent 
calculations  were  made  using  this  new  subroutine. 
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J-Factor  Calculations  Using  AGMA  908-B89 


Several  assumptions  were  made  so  that  the  equations  shown  in  AGMA  908- B89  could  be 
used  successfully  with  the  optimizing  routines.  First,  we  decided  to  use  a  rack  cutter.  There¬ 
fore,  the  number  of  teeth  on  the  cutter  was  set  at  10  000.  The  cutter  addendum  was  set  at  1.4, 
and  there  was  no  offset  for  the  cutter-tip  radius.  Second,  the  cutter-tip  radius  was  calculated  by 
using  the  equation  from  reference  3,  as  shown  in  reference  1. 

Figure  3  shows  values  of  the  J-factor  as  calculated  by  the  new  subroutine  for  the  same 
input  conditions  as  in  figure  2.  This  time,  values  remained  less  than  1  for  all  trials.  (Note  the 
scale  change.) 


Final  Runs 


The  modified  SPUROPT  program  was  run  again  with  the  same  input  conditions  as  for 
figure  1  (table  I),  and  the  process  plots  of  the  design  variables  iooked  similar  to  those  of  the 
original  run.  Final  optimal  values  for  design  variables,  noted  in  table  IV,  were  60  teeth,  2.12 
module  (12.0  diametral  pitch),  and  7.1  mm  (0.28  in.)  face  width;  fatigue  life  was  1985  million 
revolutions. 


Changing  the  shaft  diameter  from  203  to  51  mm  (8.0  to  2.0  in.)  did  not  change  optimal 
values,  although  the  tooth-bending  stress  was  lowered  slightly.  However,  when  the  program  was 
run  with  the  allowable  bending  stress  lowered  to  414  MPa  (60  000  psi)  and  the  module  limit  set 
at  4.23  (6.0  diametral  pitch)  to  permit  the  lower  limit  on  the  number  of  teeth  to  be  reached, 
optimal  values  of  the  design  variables  changed.  The  process  plot  looked  somewhat  different,  as 
shown  in  figure  4.  This  time,  the  number  of  teeth  was  driven  to  the  lower  limit,  and  the 
optimal  values  -  as  indicated  in  table  V  -  were  30  teeth,  4.23  module  (6.0  diametral  pitch),  and 
7.1  mm  (0.28  in.)  face  width,  with  a  fatigue  life  of  1564  million  revolutions.  Note  that  the  face 
width  is  a  result  of  the  limit  placed  on  facc-width-to-diameter  ratio,  and  that  center  distance 
and  bending-stress  values  were  also  at  their  constraint  limits. 


Changing  the  input  value  of  the  measured  error  on  the  pair  of  mating  teeth  affected  the 
final  optimal  values.  A  lower  measured  error  resulted  in  a  slightly  higher  number  of  teeth. 
Investigation  of  this  effect  led  to  the  discovery  of  some  small  errors  in  the  program;  however, 
they  tended  to  compensate  for  each  other.  After  these  errors  were  corrected,  the  program  was 
run  several  more  times  with  the  same  initial  gear  set  -  first  with  the  same  initial  values,  then 
with  different  bounds  and  different  objective  functions.  Run  22  used  the  corrected  program  with 
the  same  input  data  that  were  used  for  table  V;  results  are  shown  in  table  VI.  A  comparison  of 
tables  V  and  VI  indicates  that  the  program  errors  were  minor,  since  final  computed  values  were 
only  slightly  different. 

Results  for  subsequent  runs  are  given  in  table  VII,  which  shows  "standard”  bounds  for 
each  variable,  initial  values  for  the  problem,  and  resulting  t;‘im,J  values  for  design  variables. 
For  runs  in  which  the  problem  was  changed  slightly,  only  the  new,  or  changed,  bounds  are  given. 
(All  other  bounds  remained  at  standard  values.)  The  program  was  run  to  solve  three  problems: 
fatigue-life  maximization,  dynamic-load  minimWad'nn  and  i-hr  minimization.  A  special  life- 
optimization  run  was  also  made.  These  runs  are  discussed  in  the  following  sections. 
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Fatigue-life  maximization  -  Four  runs  with  different  bounds  gave  different  optima! 
results.  In  table  VII  run  22  shows  the  optimal  design  parameters  for  the  life  problem  with 
"standard”  bounds.  These  limits  (as  well  as  the  initial  values)  are  the  same  as  those  for  the 
sample  problem  in  reference  1,  except  that  the  bending  stress  is  414  MPa  (60  000  psi)  instead  of 
1241  MPa  (180  000  psi)  and  the  module,  4.23  (6.0  pitch)  instead  of  3.18  (8.0  pitch).  Fatigue  life 
increased  from  the  initial  666  million  revolutions  to  1592  million  revolutions.  Constraints  that 
were  at  their  limiting  value  are  bending  stress,  center  distance,  and  face-width-to-diameter  ratio. 
This  run  produced  virtually  the  same  design  as  noted  previously  -  before  the  small  corrections  to 
the  program  -  31  teeth,  4.03  module  (6.3  pitch),  and  7.1  mm  (0.28  in.)  face  width. 

In  run  23,  the  bounds  on  the  width-to-diameter  ratio  were  changed  to  allow  the  face  width 
to  reach  the  full  limit  of  10.2  mm  (0.400  in.).  This  resulted  in  a  fatigue  life  of  7096  million 
revolutions.  Bending  stress  and  center  distance  remained  at  their  limiting  values. 

For  run  25,  the  face-width  ratio  was  raised  to  further  observe  the  effect  of  this  parameter, 
and  actual  face-width  limits  were  changed  to  match.  This  change  resulted  in  a  very  large 
increase  in  fatigue  life,  154  475  million  revolutions.  For  this  optimization,  constraints  on  weight, 
face-width  ratio,  and  dimetral  pitch  were  at  their  limits. 

Finally,  in  run  26,  the  face-width  ratio  limit  was  increased  to  allow  a  face  width  of  up  to 
50  percent  of  the  gear  diameter.  Resulting  optimal  design  variables  were  49.2  teeth,  2.12  module 
(12  pitch),  and  a  30.7-mm  (1.21  in.)  face  width  (ratio  only  0.29);  the  fatigue  life  was  over 
200  000  million  revolutions.  Constraints  on  weight  and  pitch  were  at  their  limits. 

The  program  was  then  run  for  life  maximization  to  determine  the  effect  of  using  different 
input  values  for  the  measured  error.  Run  22  shows  the  results  of  using  a  typical  error  value  of 
0.025  mm  (0.001  in.),  as  given  in  table  I;  run  20  used  a  smaller  error  of  0.012  mm  (0.0005  in.), 
and  run  21,  a  larger  error  of  0.05  mm  (0.002  in.).  The  smaller  error  resulted  in  a  higher  number 
of  teeth  being  optimal  for  the  design  (40,  versus  31  wi<:h  the  typical  error),  with  the  same 
constraints  at  their  limits  and  about  the  same  fatigue  life.  With  the  larger  error  (run  21), 
however,  the  optimizer  was  unable  to  find  a  viable  design  -  that  is,  a  set  of  design  variables  that 
would  satisfy  all  constraints.  The  bending  stress  could  not  be  held  to  less  than  the  upper  limit, 
even  with  a  low  number  of  teeth  (minimum  diametral  pitch),  a  maximum  face-width  ratio,  and 
a  maximum  center  distance.  So,  for  run  24.  the  allowable  face-width  ratio  was  changed  to 
permit  a  full  face  width  of  10.2  mm  (0.400  in.).  This  resulted  in  a  solution  in  which  the  number 
of  teeth,  module,  center  distance,  and  face  width  were  all  at  their  limiting  values,  whereas  the 
bending  stress  was  lowered  slightly,  to  only  8  percent  over  limit. 

To  find  a  viable  solution  for  the  problem,  changes  were  made  to  run  35.  The  maximum 
face-width-to-diameter  ratio  was  changed  from  0.10  to  0.20,  and  the  maximum  tooth  width  was 
raised  to  63.5  mm  (2.5  in.).  This  resulted  in  a  solution  with  the  width-to-diameter  ratio,  gear 
weight,  and  bending  stress  at  their  limiting  values.  Optimal  design  variables  were  36.3  teeth, 
3.26  module  (7.8  pitch),  and  23.7  mm  (0.93  in.)  face  width,  with  a  fatigue  life  of  about  3000 
million  revolutions.  Thus,  with  the  higher  measured  error,  wider  teeth  were  required  to  keep  the 
bending  stress  within  limits  while  satisfying  all  other  constraints. 

Dy-fn n  ininL.ii/.quon.  -  Using  the  same  sUmdaru  bounds  as  for  <.u«,  fatigue-in'e  maxi¬ 

mization  studies,  the  SPUROPT  program  was  run  again  to  seek  the  optimal  values  of  design 
variables  with  which  to  obtain  a  minimum  dynamic  load  for  the  initial  gear  set  of  table  I. 

Results  are  noted  for  run  27,  where  the  pitch  and  number  of  teeth  were  at  minimum  limits,  and 
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the  bending  stress  and  center  distance  were  at  maximum  limits.  These  results  are  close  to  the 
optimal  values  determined  for  the  life-maximization  problem,  a  fact  not  really  surprising, 
because  of  the  relationship  between  load  and  gear  life. 

Weight  minimization.  -  Again,  with  the  same  standard  bounds  as  for  previous  studies, 
SPUROPT  was  used  to  find  the  optimum  design  values  with  which  to  obtain  minimum  weight 
for  the  initial  gear  set.  Results  given  for  run  28  show  the  number  of  teeth,  the  bending  stress 
and  the  face-width  ratio  were  at  their  limits.  To  determine  if  another  solution  with  wider  teeth 
existed,  run  30  was  made  with  a  higher  ratio  limit.  Such  a  solution  was  not  found,  however,  and 
the  minimum  weight  remained  at  0.66  kg  (1.45  lb).  It  is  not  surprising  that  these  results  are 
close  to  the  dynamic  load  optimal  values,  since  dynamic  load  is  influenced  by  the  effective 
masses  of  the  pinion  and  gear  blanks,  as  noted  in  reference  1. 

Special  life-optimization  run.  -  A  small  modification,  allowing  transmitted  power  to 
remain  constant  for  all  designs  during  optimization,  was  made  to  the  program,  so  that  we  could 
determine  if  optimal  values  would  be  affected.  (All  previous  calculations  assumed  a  constant 
force  transmission  across  all  designs.)  For  the  initial  gear  set,  the  modification  did  not  greatly 
influence  final  design.  The  optimal  results  for  run  33  were  virtually  the  same  as  for  run  25,  with 
the  same  bounds.  The  51.5  kW  (69  hp)  input  for  run  33  was  the  power  associated  with  table  I 
input  data.  Run  32  (table  VII)  shows  the  optimal  results  when  the  limit  on  the  weight  was 
raised  to  4.54  kg  (10.0  lb).  In  this  run,  the  number  of  teeth,  pitch,  center  distance,  and  weight 
were  all  at  their  maximum  limits.  Bending  stress  was  under  its  limit,  and  fatigue  life  was 
essentially  infinite  (over  10  000  yr). 


CONCLUDING  REMARKS 

Even  though  the  optimizing  routines  were  able  to  produce  viable  designs  and  provide 
optimal  values  for  design  variables,  these  values  should  not  be  considered  final  solutions  (i.e., 
not  final  design  values).  The  analytical  model  is  still  somewhat  simplified,  so  many  designs  can 
be  checked.  Each  optimal  solution  will  have  to  be  examined  more  thoroughly  for  a  given 
application. 

The  Modified  Feasible  Directions  method,  using  polynominal  interpolation  with  bounds  for 
the  one-dimensional  search,  was  the  optimization  technique  used  for  these  runs;  this  is  the 
default  algorithm.  Other  choices  were  available  in  the  program  but  were  not  used  in  this  study. 

Note  also  that  since  SPUROPT  considers  all  variables  to  be  continuous,  often  the  optimal 
value  for  the  number  of  teeth  will  not  be  a  whole  number.  The  nearest  whole  number  could  be 
chosen,  however,  and  the  program  run  again  without  using  the  optimizer  to  obtain  final  values 
for  the  calculated  parameters. 


SUMMARY  OF  RESULTS 

Computer  program  SPUROPT  was  run  for  a  set  of  spur  gear  conditions  to  determine  an 
optima!  design,  and  in  this  manner,  evaluate  the  program.  The  original  program  was  updated 
by  installing  a  subroutine  that  calculated  the  J-factor  based  on  AGMA  908-B89.  With  this 
change,  the  following  results  were  obtained: 
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1.  The  SPUROPT  program  produced  viable  designs  (i.e.,  the  objective  function  was 
maximized,  or  minimized,  while  other  parameters  were  held  within  the  bounds  of  stated 
constraints,  unless  such  a  solution  was  not  possible). 

2.  Changing  the  shaft  size  from  203  to  51  mm  (8.0  to  2.0  in.)  did  not  affect  final  results. 

3.  Changing  the  transmission  error  from  0.025  to  0.013  mm  (0.001  to  0.005  in.)  produced 
a  design  with  a  lower  module  (higher  diametral  pitch)  and  a  higher  number  of  teeth. 
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TABLE  I.  -  TYPICAL  GEAR  INPUT  DATA  FOR  INITIAL  DESIGN 

_ GEAR  AND  PINIONa _ 

Input  parameter 


Number  of  teeth  .  36 

Module,  mm  (diametral  pitch,  in.'1)  .  3.18  (8.00) 

Face  width,  mm  (in.)  .  6.4  (0.250) 

Pressure  angle,  deg  .  20.0 

Normal  load,  N  (lb)  .  1718  (386.3) 

Pinion  speed,  rpm .  5000 

Density,  kg/mm3  (lb-sec2/in.4)  .  7.835xl0"6  (7.23xl0‘3) 

Young’s  modulus,  GPa  (psi)  .  209  (30xl06) 

Poisson’s  ratio  .  0.300 

Surface  roughness,  /jm  .  0.41  (16) 

Measured  error,  mm  (in.)  .  0.025  (0.001) 

Shaft  length,  mm  (in.)  .  508  (20.0) 

Shaft  diameter,  mm  (in.)  .  50.8  (2.0) 

Brinell  hardness  . 620 

Life  factor  .  1.0 

Hardness  ratio  .  1.0 

Reliability  factor  .  1.25 

Pressure-vLcosity  coefficient,  mmZ/N  (in.2/lb)  .  0.0189  (0.0013) 

Viscosity  at  ambient  temperature,  cP  (lb-sec/in.2)  .  124  (0.18xl0‘5) 

Temperature- viscosity  coefficient,  K  (°R)  .  3900  (7100) 

Ambient  temperature  K  (°F)  .  322  (120) 


Engineering  units  are  shown  in  parentheses. 
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TABLE  II.  -  INPUT  DATA  SET  FOR  SPUROPT  PROGRAM 


$  BLOCK  A: 

OPTIMIZATION  DESIGN  FOR  LIFE  OF  SPUR  GEARS 
$  BLOCK  B: 

2.3 

S  BLOCK  C: 

0,5,7, 1000 
$  BLOCK  D: 

$  BLOCK  E: 

$  BLOCK  F: 

5 , 1 , +1 . 0 
$  BLOCK  G: 

30. ,60. ,36. 

6  .  , 12 .  , 8  . 

0.2,0.4,0.25 
$  BLOCK  H: 

1.3 

1.4 

2.5 

3.6 

3.7 

$  BLOcK  I : 

10 

$  BLOCK  J: 

9.10.1 

0 . , 0 , 60000 . , 0 

11.11.1 

0. , 0, 1.E15, 0 

12,12,1 

-l.E+10,0,0.  ,0 

13.13.1 

0 .  ,0,275.  ,0 

14.14.1 

1 .  , 0 , 10  .  ,  0 

15.15.1 

1 .  ,  0 , 3  .  ,  0 

16.17.1 

1 .  ,  0, 4 . 5, 0 

18.18.1 

3 . 5 , 0 , 5  .  ,  0 

19.20.1 

0.044,0,0.056,0 

21.21.1 

0. , 0, 1.E15, 0 
$  BLOCK  R: 

END 

36. .  36. .8. .0.25.0.25.20. 

386 . 3 , 5000 . 

0.000733,0.000733,30000000. ,30000000., 0.3, 0.3, 16., 16. 
0.001,20. ,  20. , 2 .  , 2  . 

620. .  1. . 1. . 1. . 1.25 
0.00013, 0.0000018, 7100., 120. 
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TABLE  III.  -  SPUR  GEAR  LIFE-OPTIMIZATION  RESULTS: 
ORIGINAL  PROGRAM  USING  GEOFAC 


(a)  SI  units 


Parameters 

Initial  value 

Bounds 

Final  value 

Number  of  teeth 

36 

30  to  60 

58.8 

Module,  mm 

3.18 

3.18  to  2.12 

2.15 

Face  width,  mm 

6.14 

5.1  to  10.2 

7.1 

Width- to-diameter  ratio 

0.056 

0.044  to  0.056 

a0.056 

Weight,  kg 

0.51 

0.45  to  2.01 

0.70 

Center  distance,  mm 

114 

90  to  127 

a127 

Dynamic  load,  N 

5075 

5293 

Bending  stress,  MPa 

910 

<1241 

0.004 

Contact  stress,  MPa 

1006 

<1310 

896 

Fatigue  life,  Mrev 

666 

1962 

(b)  Engineering  units 


— 

Parameters 

Initial  value 

Bounds 

Final  value 

Number  of  teeth 

36 

30  to  60 

58.8 

Diametral  pitch,  in.'1 

8.0 

8.0  to  12.0 

11.8 

Face  width,  in. 

0.25 

0.2  to  0.4 

0.28 

Width-to-diameter  ratio 

0.056 

0.044  to  0.056 

a0.056 

Weight,  lb 

1.13 

1.0  to  4.5 

1.55 

Center  distance,  in. 

4.5 

3.5  to  5.0 

a5.0 

Dynamic  load,  lb 

1141 

1190 

Bending  stress,  ksi 

132 

0.6 

Contact  stress,  ksi 

146 

130 

Fatigue  life,  Mrev 

666 

1962 

Parameter  at  constraint  limit. 
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TABLE  IV.  -  SPUR  GEAR  LIFE-OPTIMIZATION  RESULTS 

USING  JFACTR 


(a)  SI  units  with  1241  MPa  bending-stress  limit 


Parameter 

Initial  value 

Bounds 

Final  value 

Number  of  teeth 

30 

30  to  60 

a60.0 

Module,  mm 

3.18 

3.18  to  2.12 

a2.12 

Face  width,  mm 

G.  14 

5.1  to  10.2 

7.1 

Width-t  o-diameter  ratio 

0.056 

0.044  to  0.056 

a0.056 

Weight.,  kg 

0.5 1 

0.45  to  2.01 

0.70 

Center  distance,  mm 

1 14 

90  to  127 

a127 

Dynamic  load,  N 

5075 

5293 

Bending  stress,  MPa 

614 

758 

Contact  stress,  MPa 

1006 

896 

Fatigue  life,  Mrev 

666 

1985 

(b)  Engineering 

units  with  180  ksi  bending-stress  limit 

Parameter 

Initial  value 

Bounds 

Final  value 

Number  of  teeth 

36 

30  to  60 

a60.0 

Diametral  pitch,  in."1 

8.0 

8.0  to  12.0 

a12.0 

Face  width,  in. 

0.25 

0.2  to  0.4 

0.28 

Width-t  o-diameter  ratio 

O.t/56 

0.u44  to  0.056 

a0.056 

Weight,  lb 

1.13 

1.0  to  4.5 

1.55 

(.’enter  distance,  in. 

4.5 

3.5  to  5.0 

a5.0 

Dynamic  load,  lb 

till 

1 190 

Bending  stress,  ksi 

89 

<180 

no 

Contact  stress,  ksi 

146 

<  190 

130 

Fatigue  life,  Mrev 

666 

— 

1985 

''Parameter  at  constraint  limit. 
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TABLE  V.  -  SPUR  GEAR  LIFE-OPTIMIZATION  RESULTS 
USING  JFACTR  AND  LOWER  BENDING-STRESS  LIMIT 


(a)  SI  units  with  414  MPa  bending-stress  limit 


Parameter 

Initial  value 

Bounds 

Final  value 

Number  of  teeth 

36 

30  to  60 

a30.0 

Module,  mm 

3.18 

4.23  to  2.12 

a4.23 

Face  width,  mrn 

6.14 

5.1  to  10.2 

7.1 

VVidth-to-diameter  ratio 

0.056 

0.044  to  0.056 

a0.056 

Weight,  kg 

0.51 

0.45  to  2.01 

0.70 

Center  distance,  mm 

114 

90  to  127 

a  1 27 

Dynamic  load,  N 

4435 

4831 

Bending  stress,  MPa 

538 

<414 

414 

Contact  stress,  MPa 

1006 

<1310 

903 

Fatigue  life,  Mrev 

666 

1564 

(b)  Engineering  units  with  60  ksi  bending-stress  limit 


Parameter 

Initial  value 

Bound.' 

Final  value 

Number  of  teeth 

36 

30  to  60 

a30.0 

Diametral  pitch,  in.’1 

8.0 

6.0  to  12.0 

a6.0 

Face  width,  in. 

0.25 

0.2  to  0.4 

0.28 

W'idth-to-diameter  ratio 

0.056 

0.044  to  0.056 

a0.056 

Weight,  lb 

1.13 

1.0  to  4.5 

1.55 

Center  distance,  in. 

4.5 

3.5  to  5.0 

a5.0 

Dynamic  load,  lb 

997 

1086 

Bending  stress,  ksi 

78 

<60 

60 

Contact  stress,  ksi 

146 

<190 

131 

Fatigue  life,  Mrev 

666 

1564 

aPararneter  at  constraint  limit. 
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TABLE  VI.  -  SPUR  GEAR  LIFE-OPTIMIZATION  RESULTS  OF 
TAB1  E  V  WITH  ERROR  CORRECTED 


(a)  SI  units  with  414  MPa  bending-stress  limit 


Parameter 

Initial  value 

Bounds 

Final  value 

Number  of  teeth 

36 

30  to  60 

31.3 

Module,  mm 

3.18 

4.23  to  2.12 

4.03 

Face  width,  mm 

6.14 

5.1  to  10.2 

7.1 

Width-to-diameter  ratio 

0.056 

0.044  to  0.056 

a0.056 

Weight,  kg 

0.51 

0.45  to  2.01 

0.70 

Center  distance,  mm 

114 

90  to  127 

a127 

Dynamic  load,  N 

Bending  stress,  MPa 

4288 

517 

<414 

4688 

a414 

Contact  stress,  MPa 

1006 

<1310 

903 

Fatigue  life,  Mrev 

666 

1592 

(b)  Engineering  units  with  60  ksi  bending-stress  limit 


Parameter 

Initial  value 

Bounds 

Final  value 

Number  of  teeth 

36 

30  to  60 

31.3 

Diametral  pitch,  in.'1 

8.0 

6.0  to  12.0 

6.3 

Face  width,  in. 

0.25 

0.2  to  0.4 

0.28 

Width-to-diameter  ratio 

0.056 

0.044  to  0.056 

a0.056 

Weight,  lb 

1.13 

1.0  to  4.5 

1.55 

Center  distance,  in. 
Dynamic  load,  lb 

4.5 

964 

3.5  to  5.0 

a5.0 

1054 

Bending  stress,  ksi 

75 

<60 

a60 

Contact  stress,  ksi 

146 

<190 

131 

Fatigue  life,  Mrev 

666 

1592 

aParameter  at  constraint  limit. 
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Pending  Plash  < -on  tact  Mass.  Onlrr  Inter-  Width-tu-  Number  Karr  Module, *  Fatigue  Contact  Dynamic 

Hires*,  temper-  ratio  kg  distance,  ference  diameter  of  width,4  nun  lie,  Btreaa,  load 

MPa  alure,  mm  ratio  teeth4  mm  Mrev  MPa  N 
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Parameter  over  -onstraint  In 
Input  aet  at  Si  5  kW. 
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‘Deaign  variable. 

'Parameter  at  constraint  limit 
Parameter  over  constraint  limit, 
’input  «et  at  69  hp 


Face  width,  in.  Number  of  teeth 


T rial  number 


Figure  2  — J-factor  as  a  function  of  optimizer  trial 
number,  using  subroutine  GEOFAC, 


Figure  3. — J-factor  as  a  function  of  optimizer  trial 
number,  using  subroutine  JFACTR 


(a)  Number  of  teeth. 

(b)  Diametral  pitch 

2700x10® 

Trial  number 


Trial  number 
(d)  Gear  life. 


(c)  Face  width  (d)  Gear  life. 

Figure  4. — Design  variables  and  gear  life  as  a  function  of  optimizer  trial  number  for  allowable  bending  stress  of  414 
MPa  (60  000  psi). 
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